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Abstract 



Within Current Density Functional Theory, we have studied a quantum 
dot made of 210 electrons confined in a disk geometry. The ground state of 
this large dot exhibits some features as a function of the magnetic field (B) 
that can be attributed in a clear way to the formation of compressible and 
incompressible states of the system. The orbital and spin angular momenta, 
the total energy, ionization and electron chemical potentials of the ground 
state, as well as the frequencies of far-infrared edge modes are calculated as 
a function of B, and compared with available experimental and theoretical 
results. 

PACS 73.20.Dx, 72.15.Rn, 78.20.Bh 
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I. INTRODUCTION 



The study of the formation of compressible and incompressible states of quantum dots 
has attracted an enormous interest in the last few years. It has been argued that for 
a sufficiently smooth confining potential, the edge of the dot exhibits to a certain degree, 
narrow incompressible regions where the electronic density is determined by the filling factor 
v as p = v/(2irC 2 ), where C=(hc/eB) 1 ^ 2 is the magnetic length, and compressible regions 
where p may vary with the radial distance r. These states are at the origin of phenomena 
such as, e.g., "magic" angular momentum quantum numbers, composite droplets and edge 
reconstruction ||||. Some evidence for the formation of these states has been found in the 
magnetic field (B) dependence of the addition energies, i.e., the energies to add one more 
electron to the dot [0,|J, an d more recently, in the oscillations of the resonance frequencies 
of edge magnetoplasmon modes of dots and antidots . 

Exact diagonalizations of the electronic Hamiltonian corresponding to dots with a small 
number of electrons N (N < 8) [P~0| [T~3|] have shown that these systems have strongly corre- 
lated ground states (g.s.), especially at high magnetic fields. Approximation schemes such 



as the Constant Interaction model of Ref. 141, the Self- Consistent model of Ref. 15 , the 



Hartree-Fock method of Refs. |]6|JT6[1, or other mean field approximations [17] which have been 
used to carry out calculations in larger dots, do not take into account these correlations, so 
they may fail in reproducing a number of g.s. features that appear as a function of B ||T3 



It is worth to mention that all these approaches have been only applied to the description of 
medium size (N <100) dots and relatively small values of B due to the technical difficulties 
one has to face otherwise because of the high density of single particle (s.p.) states. 

Recently, Current Density Functional Theory (CDFT) [|18f has been used to study g.s. 



properties of small quantum dots [19]. The accuracy of the method has been tested in the 



case of two and three electron systems, for which exact calculations are available It 



has been shown [19| that CDFT reproduces accurately the exact results, thus providing a 



convincing evidence that it is a powerful tool to quantitatively describe the g.s. properties 
of quantum dots in magnetic fields. CDFT has also been used to study the far-infrared edge 
magnetoplasmon modes in medium size quantum dots by means of an equation of motion 
method J2D]. 

In the present work, we extend the CDFT calculations to the case of a large dot of radius 
i?=1600A containing 210 electrons in a magnetic field ranging from B=0 to ~ 15 T. The 



far-infrared spectrum of this dot has been experimentally studied in [21]. The large number 



of electrons makes this dot closer to the 2D conditions than any other system previously 
considered in a calculation of this kind. Moreover, the high 5's here attained allow us to 
study g.s's. which are associated in a clear way with bulk incompressible states with filling 
factors v — |, and from u—1 to 10, and to several incompressible states in between. Such a 
rich structure has not been previously disclosed in a single large dot, where genuine magnetic 
effects, even at rather small B, are not masked by surface nor finite size effects. 

The density functional we are using does not incorporate one of the key features to 
describe the occurence in a single dot of several incompressible strip regions in the fractionaly 
quantum Hall regime, namely cusps in the exchange-correlation scalar potential at certain 
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filling factors (see for example Refs. 0,0). However, it is able to describe the appearance 
of bulk regions with a fractionary filling factor. This is illustrated in the case of the v = | 
incompressible state. 

To fully display these magnetic field effects, we have performed the calculations at a 
very small temperature (T), 0.1 K. The use of such temperatures in the calculations is 
a challenge for large dots and intense magnetic fields. Some results at higher T are also 
discussed. Finally, we have obtained the far-infrared spectrum of the dot. 

II. CURRENT DENSITY FUNCTIONAL DESCRIPTION OF QUANTUM DOTS 



Current Density Functional Theory is comprehensively described in ||18|| , and has been 
already applied to quantum dots [|i9p0f . In the present work, we shall be considering N 
electrons moving in the z=0 plane where they are confined by the dot potential V+(r) with 
r = y/x' 2 + y 2 . On this system it may act a constant magnetic field in the z direction 
described by the vector potential A = ^(—y,x,0)B. Introducing the cyclotron frequency 
ijj c = eB/mc, we can write the CDFT total grand potential A = E — TS — /iiV adding to 
the total energy 

B =i/dFr(f) + f jifrm + U /*V#) +S >»BE/, Ss 

i 

(1) 

+ J drV + {r)p{r) + l - J J dfdf' + E xc (p, £, V) , 

the temperature times the total entropy S given by: 

S=-£[/iln/i + (l-/*)h(l-/<)] . (2) 

i 

and the electronic chemical potential \i times the electron number N = J2ifi- 

To easy the formulae, we have used effective atomic units, defined by h = e 2 /e = m = 
kg =1, where is the Boltzmann's constant and m = m*m e is the electron effective mass. 
In this system of units, the length unit is the Bohr radius times e/m*, and the energy 
unit is the Hartree times m*/e 2 , which we call respectively, ag and E* H . m* is the electron 
effective mass in units of the bare one, which together with a dielectric constant e and a 
gyromagnetic factor g* are characteristics of the semiconductor (for example, g* = —0.44, 
e=12.4 and m*=0.067 in GaAs), and fiB = he/(2m e c) is the Bohr magneton |[23|| . For GaAs 
we have a* ~ 97.94 A and E* H ~11.86 meV. 

In Eqs. (0,0), s Zi is the z-component of the spin, and fi is the occupation number of the 
i-s.p. level. The particle p(r ), kinetic energy r(r ) and paramagnetic current densities j p (r) 
are defined in terms of the occupation numbers and s.p. wave functions 0j as follows: 

P(r) = £/<IM0I 2 (3) 

i 

r(r)=E^|V0,(r)| 2 (4) 
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--EM 



(5) 



where eg is the azimuthal unit vector. Due to the axial symmetry, the electronic wave 
functions can be factorized as 



e- M u nla (r) 



(6) 



where a = ±1, n— 1,2,3..., and I = 0, ±1, ±2, ±3... These wave functions are eigenstates 
of the angular momentum around the z-axis with eigenvalue —I. At B ^ 0, the s.p. level 
i is non degenerate, and in the above expressions we have used the short notation i = 
{n, I, a}. V + (r) is the Coulomb potential generated by a jellium disk of radius R and areal 
density p + = p&6(i? — r) p4| , and the ph value has been fixed by charge neutrality, i.e., 
npbR 2 = N+ = N. For the i?=1600A and iV=210 system which corresponds to one of the 
GaAs dots studied in ETH, one has p b =0.239 a* ~ 2 . The use of a step function to model 
the positive neutralizing background, i.e., a sharp confining potential fl25| , together with a 



smooth exchange-correlation potential, implies that no strips of filling factor smaller than 
the one at the center of the droplet are expected to appear at the edge of the dot. That will 
be confirmed by the calculations presented in the next section. 

The exchange-correlation energy E xc (p,£,V) in Eq. ([l|) has been taken from [19]. It is a 
functional of the particle density p, the local spin polarization £ 



pt - pi 

PT + Pl 



(7) 



where p^(pi) is the particle density of spin up (down) electrons, and of the local vorticity 
V(r) 



V(r) = V(r) 




(8) 



For further reference, we also define the local magnetization m(r) as m(r) = p(r) £(r). 

To obtain the s.p. wave functions 4>i{r) and occupation numbers fi, one has to minimize 
A = A —J2i^i(4 ) i I 4>i)- Minimizing A with respect to the wave functions one obtains the 
Kohn-Sham (KS) equations: 



1 / d 2 Id 

2 \ dr 2 r dr 



l + -uy + V+(r) 



(9) 



+ / dr 



/ p(r') 



e , A xc [r) . . 1 

-I h V xca {r) + -g p B Bo 

c r I 



with £j = Xi/ fi and where 



V xca (r) 



SE XC (P,Z,V) 



Spa 



^ -A xc (r) ■ k 



(10) 
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being A xc the exchange-correlation vector potential 

e t i \ e a - c d ( &E XC \ 

-A^ r = -A^ea = — — — — — L f e e . 11) 

c w c epdr V 5V lp '7 v ; 

Minimizing A with respect to and making use of the KS equations, one obtains the 
occupation numbers: 

f* = 1 + e (ei-»)/T ■ ( 12 ) 

The electron chemical potential is fixed from the normalization condition: 

"-S l + elU/r - < 13 > 

The KS differential equations and the normalization condition Eq. (|T3|) have been solved 
selfconsistently, without expanding the s.p. wave functions in a necessarily truncated basis 
of Fock-Darwin or Landau orbitals. Our iterative method works for weak and strong B fields 
as well, for which the effective potential is very different. It has the advantage of avoiding 
the study of how the results depend on the size of the basis. 



III. GROUND STATE RESULTS 

We have studied the g.s. of the iV=210 electrons dot at T=0.1 K varying the magnetic 
field from to ~15.5 T in order to explore states with filling factors v = | and above. We 
recall that integer z/s also represent the number of occupied Landau bands, which we label 
as (M, |) or (M, j), depending on the spin of the s.p. states in the band. M = n + (\l\ — l)/2 
is the Landau level index. 

The density and single particle energies corresponding to the B=0 case are represented 
in Fig. ([I]). The s.p. energies are distributed along parabolic-like curves as a function of I, 
each curve being characterized by a different value of n. The figure shows the well known 
fourfold degeneracy of each single particle level, corresponding to the possibilities ±1, (T-l)- 
It is remarkable that even for such a large system, the density still displays many oscillations 
in the bulk which have to be gradually washed out as the dot evolves towards the 2D system. 

Figure (§ shows the CDFT s.p. energy spectra for B= 10.28, 5.14, 3.43, 2.57, 1.71, and 
1.29 T as a function of n and Observe that due to the sign of the linear /-term in Eq. @, 
most occupied levels have now a positive I value, i.e., a negative orbital angular momentum, 
and that the (M, |) band lies below the (M, 1) one because of the negative value of g*. 
Fig. (|3p shows the corresponding g.s. particle p(r) and local magnetization m(r) densities 
together with the predicted values of total orbital and spin angular momentum (L z , 2S Z ). 

From these figures, one can see that these states are the finite-size analog of the 2D 
Landau incompressible states with v=\ to 4, 6 and 8 which, for an electronic density p e = 
P6=0.239 clq~ 2 , would be precisely realized at the above B values. Indeed, the densities in 
Fig. (|3|) exhibit a step-like shape whose plateau has a density determined by an integer 
filling factor v as p = veB /2irc. To help follow the discussion of the results we present, we 
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have collected in Table I the values of B and C that correspond to v=\ to 10. Notice that 
even for i/=10, R is much larger than C, so that the filling factor v can be safely used to 
describe the occupation of the Landau bands in the interior of the droplet. 
The g.s. values of the orbital angular momentum agree well with the law 

L = ^ NiN-J) (M) 
and the spin is approximately given by 

expressions that are valid in the big-iV limit for the Landau bands. This is explicitly shown 
in Fig. (|4]), where we have plotted the calculated orbital and spin angular momenta as a 
function of B. The small plateaus along the L Z (B) curve at integer v, also visible in the 
2S Z (B) curve, show the robustness of these states against changes in B, and their tendency 
to not to change their (L z , S z ) content. 

In panel (a) of Fig. (|3|) we have plotted the density of the maximum density droplet 
(MDD) state |f|: 

1 N ~ x I ( r \ 2 ' _ . 2 
PMnn(r) = — £ - e~** . (16) 

One can see that both densities are very similar, the only difference appears in the edge 
region. The lump in the CDFT density at the edge is due to the rather sharp confining 
potential, on the one hand, and on the other hand to the exchange-correlation energy, which 
is changing rapidly in the surface region of the droplet, see Eqs. ( |10| , |lTD , and Refs. [ p2fl . 

It is also interesting to look at the energy difference between the (M, f) and (M, [) 
Landau bands shown for example in Fig. (0). It is very small for even u, which correspond 
to g.s.'s having S z ~0, because the Zeeman term in Eq. ([]]) is also small (0.127 meV at B=5 
T, for instance). On the contrary, that separation is sizeable for the v—1 and 3 cases, in the 
2-4 meV range. These g.s. states have S z ^0. One should have in mind that u c , which gives 
the difference of frequencies between Landau levels when spin and interaction effects are 
neglected, is around 18 meV for B=10 T. V xca is the origin of the large spin splittings. The 
importance of the exchange energy in the spin splitting of the Landau bands was already 
pointed out in ||16|| . The similarity of the g.s. corresponding to v=\ and the MDD state is 
due to V xca , as it prevents the (1, |) and (1, |) Landau bands from being close at 5=10.28 
T and thus, to contribute both to the building of the g.s. 

We have not attempted to disentangle which part of the splitting we have found is due 
to the exchange energy, and which part is due to the electronic correlations. It would 
be possible at low B as the exchange-correlation energy comes from the Tanatar-Ceperly 
parametrization [f27| , but not at high B because of the form of the parametrization given by 
Levesque et al ||2£j . Exchange-correlation effects enhanced by the finite size of the system 



are also responsible for the differences between the calculated values of L z , and especially of 
S z , and the ones given by Eqs. (|14],[nj). In the case of S z , these differences can be sizeable 
for 'large' odd filling factors, see Fig. (|). 
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Figures @ and (§) show the s.p. energies and densities for B= 1.9, 3, 4, 7, 9, and 11 T. 
These values originate compressible g.s.'s in the dot, and a careful look at these figures sheds 
light about the appearance of compressible and incompressible states in the finite system. 

When the magnetic field is increased from that corresponding to v=l, some electrons are 
promoted within the (1,T) band to s.p. states with higher / values. This is seen in panel (f) 
of Fig. (H). Indeed, it can be realized from that figure that the chemical potential intersects 
the s.p. curve at a value of I larger than 210, and that at the same time, some low-/ s.p. 
levels are at an energy slightly above fi. This is precisely the evolution of the MDD state 
when B is increased and edge reconstruction occurs 



A different behavior is shown in the other panels of Fig. (^). Consider for instance, the 
situation displayed in panel (d). Around 140 electrons are in the (1,T) band, as may be 
inferred from the intersection point of that band with the \i line. All other electrons are 
oddly distributed in the (1, j) band, starting from the s.p. levels with higher orbital angular 
momentum, whose energies are well below the chemical potential line because of the dip 
around /=120 caused by the exchange-correlation energy term. Whenever the partially filled 
Landau band presents a flat region, the electron chemical potential has to be the average 
s.p. energy of that region. Otherwise, electrons will continue to fill s.p. levels with high-/ 
values belonging to the previous band. 

The corresponding electron density manifests the peculiarities of the s.p. spectrum in 
quite a distinct manner: the oscillations displayed by the density at B=7 T as shown in 
panel (d) of Fig. fl5|), are washed out at B=9 T as shown in panel (e) of the same figure 
because for this magnetic field, the s.p. energies of the (1,|) band are distributed more 
evenly around \i. 

The evolution of the compressible region between v=2 and 1 when B increases from 5.14 
to 10.28 T deserves further discussion, as this region extends in a wide B range and it makes 
the magnetic effects easier to visualize (a similar evolution has been found between u=3 and 
2). After a regime extending up to > 7 T in which the droplet densities present an oscillatory 
behavior, at B ~ 8 T the (1,|) flats out and accordingly, the density oscillations disappear. 
The configuration at B=8 T having L z = 17 031 and 2S Z = 124 cannot be identified with 
the u—1 incompressible state, since L z and 2S Z are far from having the appropriated values, 
neither other similar g.s. configurations corresponding to magnetic fields as high as 10 tesla, 
see Fig. (^j). It means that looking only at the densities might be misleading for identifying 
integer filling factor states. The evolution of the 5=10 T state towards the v=l one at 10.28 
T proceeds through an interesting change in the s.p. spectra of the (1,|) band: it bends 
upwards crossing the \i line, so that its occupied higher-/ states get now empty. To keep 
constant the number of electrons on the one hand, and to keep L z increasing with B on the 
other hand, the (l,f) band extends its occupied s.p. states up to higher-/ states. A further 
increase in B 'decouples' /i from the (1,1) band, and the situation is eventually that shown 
in panel (a) of Fig. (g). 

It is of some importance to figure out what is the quantitative influence of temperature on 
the above features of the g.s. To this end, we have carried out calculations using a moderate 
T, namely 2 K. For quantum dots as the one studied here, it appears that the density of 
states around the chemical potential is large, and such a small temperature already produces 
changes in the s.p. spectrum and the electron density as well. This is illustrated in Fig. (^), 
where we have drawn the s.p. energies and electron densities corresponding to B=4 and 7 
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T. 

A comparison of cold and hot s.p. energies in Fig. (0) shows that a first thermal effect 
is that the small wiggles around /i in the cold spectrum are progressively wiped out when 
T increases. In other words, the low-/ states which are unoccupied at T=0.1 K (they lay 
slightly above the chemical potential), are now occupied with some probability. This flattens 
the electron density in the bulk of the droplet. A second thermal effect is the deepening of 
the high-Z s.p. levels as T increases. This is due to the widening of the surface region, a 
situation that has been already met in other finite fermionic systems, see for example Refs. 
?T|-P^[. This thermal effect reduces the amplitude of the density oscillation at the edge of 



the dot. 

It is thus quite apparent that to study the incompressible states, one has to use temper- 
atures in the right range of values. Otherwise, the thermal average implicit in the definition 
of one-body densities will wash out the quantum oscillations of compressible states, making 
it more difficult the identification of the incompressible ones. 

At B= 15.42 T, an incompressible state appears, which deserves a separate analysis. 
From the relationship p = veBjl'Kc, it would correspond to a filling factor v — |. We 
predict for this state an angular momentum L z =32 511, which is very close to the one 
given by Eq. flUp, namely 32 917. We recall that the angular momentum of a droplet in 
a Laughlin's state |34j] is given by that equation only when v~ l equals 1, 3, 5,... We argue 



that this formula is a good approximation not only for integer filling factors and those of 
Laughlin's states, but also for v = |. As a matter of fact, Fig. (f|) shows that Eq. (|14"|) 
works well even for compressible states. 

As we show in panel (c) of Fig. (P), the electron density is constant in the bulk and has 
a bump at the dot edge. This has also been found in simulation calculations of Laughlin's 
states However, these simulations have been carried out for few electron dots, and 



density oscillations are still visible in the bulk. In the CDFT calculation here presented, the 
density is constant in the bulk which, together with the associated values of L z and 2S Z , 
makes it undubious to label this dot state with a filling factor v — |. 

Further insight into this state can be obtained splitting its density as follows: 

AT-l 

P\( r ) = £ fi I <M r ) I 2 + fi I <M r ) ! 2 = Pbuik(r) + psurf(r) , (17) 

i=0 i>N 

with i = {1,1, T}. This produces the result displayed in panel (c) of Fig. fl8|), where Pbuik( r ) 
is represented by a dashed line, and p sur f{r) by a dot-dashed line. The integral over r of 
the former gives N bu i k ~140, and that of the later, N sur f ~70. This means that two-thirds 
of the electron number in the dot is in the bulk region, and one third in the edge region. 
A look at panel (a) of Fig. (Q) shows that the occupation number of the s.p. levels with 
/ < 210 is quite similar, around 2/3. Of course, fi are thermal occupation numbers, and to 
have a fractional value, a finite but small T has to be used in the calculation (0.05 K in the 
present case). However, it is a workable way CDFT has to accomodate noninteger filling 
factors in the bulk of an extended droplet. Whether the occupation numbers represent a 
'simple' thermal situation or come from a genuine correlation effect, depends on whether 
correlation effects have been smeared out by an unadequate high temperature or not. We 
want to stress that the temperature we have used in the present calculations (50 or 100 
mK) is much lower than the one at which most experiments on quantum dots have been 
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carried out. It is worth to recall that within the zero temperature KS scheme, fractionary 
occupation numbers may appear for s.p. i-states such that £j = \i [p7 |. The generalization of 



Density Functional Theory to this situation is called Ensemble Density Functional Theory, 
and has been used in |22| to describe quantum dots in the fractionary quantum Hall effect 
regime. 

It is worth to point out that confining potentials produced by a jellium disk favor the 
existence of incompressible regions in the bulk of the dot, since there is a tendency in the 
electron bulk density to have a value around that of the jellium background in order to 
screen the Coulomb potential. This can be seen for example in the compressible states 
shown in Fig. (||). It is also clear that whenever an incompressible state appears, its bulk 
density is that of the neutralizing jellium, see Figs. (§,§). However, a physically sound 
exchange-correlation energy is at the very origin of the incompressible states we have found. 

To substantiate our case, we have displayed in panels (b) and (d) of Fig. (||) the results 
obtained at 5=15.42 T when the exchange-correlation energy of is dropped from the 
functional. This supports that the results shown in panels (a) and (c) of that figure are not 
only the trivial consequence of screening the jellium density, nor of using an unadequate high 
temperature, but a product of the high-I? correlations included in the density functional. 
Actually, one can see from panel (d) that density oscillations are clearly visible even at T=2 
K. 

For integer v values, the changes are not so dramatic, and the Tanatar-Ceperly exchange- 
correlation energy alone yields results quite similar to those we have shown for v > 1. 

The B dependence of the total energy E is shown in the top panel of Fig. (|9|). Local 
minima in E are clearly visible at or near even v values corresponding to paramagnetic 
states of the droplet having 'zero' total spin momentum. Local maxima are at or near odd 
v values which correspond to ferromagnetic states having large 2S Z values, but smaller than 
these given by Eq. (0). The energy also presents an inflection point at the ferromagnetic 
v=\ state with 2S I 2 =210. In contradistinction with the situation for small dots where only 
one transition from paramagnetic to ferromagnetic g.s.'s is observed [38], several smooth 
transitions take place below B ~ 10 T for iV=210. Also shown in that figure (bottom 
panel), is the dependence of E with L z , with local minima at the L z values corresponding 
to even filling factors. 

The existence of several paramagnetic states plays a crucial role in keeping the energy of 
the droplet at around the value corresponding to B=0. \E\ decreases as S z increases, and 
viceversa. It is only from the z/=2 state on that E increases monotonously, as S z does up to 
reaching the full polarized value around v—\. 

The electron chemical potential and the ionization potential IP = E^-i — E^ 



arc 



shown in Fig. ([T0|). They display a "sawtooth" behavior with B, with a fast falling near 
integer v values. The corresponding B values or equivalently, the corresponding L z values, 
are a kind of "magic numbers" for which the dot is particularly stable. It is worth to notice 
that the large size of the studied dot is the reason why -IP and fi are so similar, and why 
the oscillations in IP are rather wide, extending over several tesla and giving them the 
possibility of being experimentally observed more easily than in small dots |7||§. 

We have not attempted a systematic study of the addition energies E a d, but have obtained 
them for the B=4 and 7 T cases. To this end, we have calculated £'(209) keeping constant pb 
and adjusting R so as to have iV + =209. Subtracting £(209) from £(210) it yields E ad ~192 
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meV and 190 meV, respectively. These values agree well with those obtained with the 
schematic model of ||13|| , which are 215 and 213 meV, respectively. 



IV. FAR-INFRARED EDGE MODES 

In this Section we present the results obtained for the far-infrared edge modes of the 



iV=210 dot and compare them with the experimental results of pl| . The method we have 
used is thoroughly described in [pO|j . It is based on a variational solution of the equation of 
motion obeyed by the excitation operator that generates the modes, much along the method 
originally proposed by Feynman to describe density excitations of superfluid 4 He [j4U|J . 

The expressions to obtain the frequencies of edge modes of different multipolarity L are 
rather cumbersome, but only make use of one-body densities associated with the ground 
state of the system. In the dipole L=l case, the one for which experimental results have 
been unambiguously obtained, we worked out a simple expression for the u±i frequency |20 



w±i = + dfAV + (r)p(r) ± ^ . (18) 

It is easy to further elaborate this equation and write the dipole frequency as a function of 
the electric field £{r) created by the electrons at the edge of the dot. We obtain: 



" ±1 = ^ + ^±f. (19) 



We have drawn £ (r) in Fig. ( |TT|) for the B values corresponding to v=\, 1, and 4. As it has 
been previously obtained p] , it is sharply peaked a the dot edge. The shoulders displayed 
by £ {f) near the edge reflect the complex structure of the droplet density in that region. 

We show in Fig. (|12|) the frequencies corresponding to the L — ±1, ±2, and ±3 modes, 
together with the experimental points of ||21|| . We want to stress that there are no adjustable 
parameters in the calculation. From the figure, one sees that the agreement between the- 
oretical and experimental spectra is good. In particular, the L — ±1 branches are well 
reproduced up to a high value of B. The interference between the u>_3 and u + i which Shikin 
et al. E2] suggest it causes the level repulsion experimentally found, is at the right place and 
has the correct amplitude. In our calculation, the quadrupole branches fit some experimental 
points. However, it is not yet clear whether these points correspond to quadrupole excita- 



tions ||43|| , or are the result of a fragmentation of the dipole spectrum, as some calculations 



in small dots might indicate j44|] 



In a recent work [EJ, it has been pointed out that incompressible strips at the edge of 
quantum dots could be detected by far-infrared spectroscopy when the confining potential is 
nonparabolic. In particular, the low-frequency dipole branch should be especially sensitive 
to the shape of the dot edge and hence, exhibits oscillations as B is varied. Eq. (|T9| ) shows 
that these oscillations appear if the electric field at the edge of the dot has an oscillatory 
behavior with B, and are absent otherwise. For example, when the confining potential V + (r) 
is parabolic, Eqs. (|18| , [T9|) yield nonoscillatory dipole frequencies, as £/R is a constant. This 
is easily seen from Eq. ([18]). In particular, that constant is oJq = npb/R if the parabolic 
approximation of V+ in the r/R « 1 limit, V+(r) = is justified. 
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We have analyzed the negative dispersion branch of the dipole mode represented in Fig. 
(|l2l), plotting the ratio u_/uj d : as a function of B, where ui c } is given by Eq. (|i"9"D taking for 
£{R) its value at B=0. The results are shown in Fig. (|13|). One can see that the oscillations 
in uj_ are indeed filling factor related, and as indicated in Ref. @, the minima correspond 
to half filled Landau levels (odd z/), and the maxima to fully occupied Landau levels (even 
v). Filling factor u—1, not attained in U, is an exception as it constitutes in a sense a filled 
level by itself. The dipole positive dispersion branch also presents oscillations, but they are 
an order of magnitude smaller than these of the u;_ branch. 

Apart from the z/=l and 2 cases, the amplitude of the oscillations we find is smaller than 
in ||, but one should have in mind that different confining potentials are used. It is also 
worth to notice that in the present calculation, the oscillations in the dipole frequency are 
not due to the existence of integer filling factor strips at the surface of the dot, which as we 
have indicated before, do not appear here because of the sharp jellium density we are using. 
Rather, they arise due to the appearance of integer filling factor regions in the bulk of the 
dot. 



V. SUMMARY 

In this work we have used Current Density Functional Theory to describe the ground 
state and multipole spectrum of a quantum dot made of 210 electrons, confined in a disk 
geometry and submitted to weak and intense magnetic fields as well. The sharp jellium 
confining potential does not allow the system to exhibit incompressible states at the edge, 
but its large size permits these states to develop in the bulk. We have been able to identify 
such incompressible states with filling factors from one to ten, and many other compressible 
states in between. 

At 5=15.42 T, a v = | incompressible state appears with two-thirds of the electron 
number in the bulk region, and one-third in the edge region. The occupation number of the 
single particle states in the bulk region is around 2/3, and near the edge they rise to unity. 
These features have also been found in Ref. [^] when the confining potential is supplied 
by a positive jellium background, and have been interpreted there as the formation of a 
composite edge in the system. 

The finite temperature version of CDFT permits to solve a technical problem when the 
density of s.p. states is large and it is impracticable to solve the zero temperature KS 
equations with s.p. occupation numbers or 1. Interestingly, at low enough temperatures 
it permits to obtain the electron density and other closely related characteristics of highly 
correlated states that cannot be described even approximately in terms of integer occupation 
numbers, provided of course, these correlations are built in the functional. This is precisely 
the present case, and also that of |^,22]. 



We predict values of the orbital and spin angular momentum that agree well with the 
ones pertaining to Landau and Laughlin states. States with integer filling factors are rather 
robust against changes in B. In particular, the study of the total energy, ionization and 
electron chemical potentials as a function of B has shown that even S z ~ paramagnetic 
states are especially stable, corresponding to local minima of the total energy, and with 
large ionization and electron chemical potentials. Odd u, ferromagnetic states with large 
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S z values correspond to local energy maxima, apart from the particular fully spin polarized 
u=l state, and also have comparatively small ionization and electron chemical potentials. 

We have stressed the important role played by the exchange-correlation term of the 
current density functional in the quantitative description of these effects on the one hand, 
and the need to perform the calculations at small temperatures to disclose fine features of 
incompressible states, on the other hand. 

Finally, we have studied the multipole spectrum of the dot. Good agreement with the 
experimental data has been found. For the dipole mode, we have been able to associate the 
oscillations with B of the dipole frequencies around the classical value, with the oscillations 
with B of the electron electric field at the edge of the dot. These oscillations are filling factor 
related and in our case are associated to the formation of bulk incompressible disks, and in 
the case of 0] to the formation of edge incompressible strips. 
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TABLES 



TABLE I. Values of B and C corresponding to integer filling factors v=\ to 10. 



V 


B (T) 


C (og) 


1 


10.28 


0.82 


2 


5.14 


1.15 


3 


3.43 


1.41 


4 


2.57 


1.63 


5 


2.06 


1.82 


6 


1.71 


2.00 


7 


1.47 


2.16 


8 


1.29 


2.31 


9 


1.14 


2.45 


10 


1.03 


2.58 
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FIGURES 



FIG. 1. Top panel: s.p. energies as a function of orbital angular momentum I at 5=0. The 
full, upright triangles represent spin up states, and the empty, downright triangles, spin down 
states. The horizontal solid line represents the electron chemical potential. Bottom panel: The 
corresponding electron density as a function of r, both in effective atomic units. The dotted line 
represents the jellium density. 

FIG. 2. Single particle energies as a function of orbital angular momentum I for v=l, 2, 3, 
4, 6, and 8. The horizontal solid lines represent the electron chemical potential. The full, upright 
triangles represent (M, f) bands, and the empty, downright triangles, (M, j) bands. In panel (a), 
only M=l bands appear, whereas in panel (f), bands with M=l to 5 are visible. 

FIG. 3. Electron p(r) (solid lines) and local magnetization m(r) (dot-dashed lines) densities 
corresponding to the configurations displayed in Fig. (2). The dashed curve in panel (a) represents 
the density of the MDD state. All densities are in effective atomic units. Also shown are the values 
of (L Z ,2S Z ). 

FIG. 4. L z (circles, left scale) and 2S Z (triangles, right scale) as a function of 5. The values 
corresponding to integer filling factors 1 to 11 are denoted by solid symbols, and the line connecting 
the 2S Z points is to guide the eye. The dashed line represents Eq. (14). 

FIG. 5. Single particle energies as a function of orbital angular momentum I for 5=1.9, 3, 4, 
7, 9, and 11 T, respectively. The horizontal solid lines represent the electron chemical potential. 
Same notation as in Fig. (2). 

FIG. 6. Electron p(r) (solid lines) and local magnetization m(r) (dot-dashed lines) densities in 
effective atomic units corresponding to the configurations displayed in Fig. (5). Also shown are 
the values of (L z , 2S Z ). 

FIG. 7. Top panels: s.p. energies at T=0.1 K (dots), and at T=2 K (triangles). Bottom 
panels: Electron density p(r) at T=0.1 K (solid line), and at T=2 K (dashed line). The left panels 
correspond to B=4 T, and the right ones to B=7 T. The electron chemical potential is indicated 
by a horizontal line (solid, T= 0.1 K; dashed, T= 2 K). 

FIG. 8. Top panels: s.p. energies (triangles, left scale), and occupation numbers (dots, right 
scale), as a function of orbital angular momentum I for 5=15.42 T. The horizontal solid and dotted 
lines represent, respectively, the electron chemical potential, and the occupation number |. Bottom 
panels: the corresponding electron density as a function of r, both in effective atomic units. In 
panel (c), the bulk density is indicated by a dashed line, and the surface density by a dot-dashed 
line. The dot-dashed line in panel (d) represents the electron density at T=2 K. Also shown are 
the values of (L z , 2S Z ). The right panels have been obtained dropping high- 5 contributions to the 
exchange-correlation energy coming from Ref. [29]. 
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FIG. 9. Top panel: Total energy as a function of B. The values corresponding to integer 
filling factors 1 to 11 are indicated by solid symbols. The line connecting the dots is to guide the 
eye. The insert shows the low-B region, with the energy scale at the right. Bottom panel: Same 
as top panel but as a function of L z . 

FIG. 10. Electron chemical potential (circles), and —IP (triangles), as a function of B. The 
lines connecting the points are to guide the eye. The points corresponding to integer filling factors 
1 to 11 are indicated by solid symbols. 

FIG. 11. Electron electric field £(r) in effective atomic units as a function of r for filling 
factors v = |, 1, and 4. 

FIG. 12. Edge mode frequencies as a function of B. At B ~ 0, from top to botton the curves 
correspond to L= ±3, ±2, and ±1. The experimental dots are from Ref. [21]. 

FIG. 13. lo^/lo 01 as a function of l/B (lower scale) and v (upper scale). The solid symbols 
correspond to integer filling factors 1 to 11. 
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